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Measurements of transition form factors for the electro-excitation of nucleon 
resonances in the Q 2 = 5 — 14 GeV 2 region can provide one with the informa- 
tion on quark wave functions at small transverse separations. In particular a 
comparison of form factors for the states of opposite parity can give insight in 
the mechanisms of chiral symmetry breaking. I discuss perspectives of the the- 
oretical description of such reactions using a combination of lattice calculations 
and light-cone sum rules. 
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1. Introduction 

Electroproduction of nucleon resonances has long been recognized as an 
important tool in the exploration of the nucleon structure at different scales. 
There is a growing consensus that perturbative QCD (pQCD) factorization 
based on hard gluon exchange is not reached at present energies; however, 
the emergence of quarks and gluons as the adequate degrees of freedom is 
expected to happen earlier, at Q 2 ~ a few GeV 2 . Measurements of the form 
factors in this transition region are planned at Jefferson Lab 1 ' 2 using the 
CLAS12 detector. In this talk such reactions arc addressed from the theory 
perspective. I try to formulate the physics goals and explain an approach 
to the electroproduction of resonances that combines lattice calculations 
of wave functions at small transverse separations with dispersion relations 
and quark-hadron duality, known as light-cone sum rules (LCSR). 

Quantum chromodynamics (QCD) predicts 3-6 that at large momentum 
transfer the transition form factors become increasingly dominated by the 
contribution of the valence Fock state with small transverse separation be- 
tween the partons. In reality the dominance of valence configurations is 
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likely to be achieved but "shrinkage" to small transverse separations does 
not seem to occur and also power counting rules based on helicity conserva- 
tion do not work at moderate Q 2 that we are able to study experimentally. 

To be more specific, the nucleon valence state contains contributions 
with different orbital angular momentum l z ~ 0, ±f . The leading contribu- 
tion to form factors at large Q 2 comes from the l z = state: 6 

| Pt) /.=o = [ [dx][*k] ^,=o ( g. )x 
J %^Jx x x 2 x 3 

x | |u f (a;i, k 1 )u l (x 2 , k 2 )d f (x 3 , k 3 )) ~ \u t {x 1 ,k 1 )d l (x 2 , k 2 )u f (x 3 , fc 3 ))| 

where ip e * =0 (xi, ki) is the three-quark light-cone wave function that depends 
on quark momentum fractions Xi and transverse momenta kf . 

The simplification that occurs at asymptotically large Q 2 is that the 
fc-dependence of wave functions becomes irrelevant and all necessary (non- 
perturbative) information is contained in the integral over transverse mo- 
menta 

*3{Xi-,n) = j\d 2 k] ^'=°(Xi,ki) 

where the cutoff ji ~ Q has to be imposed to make the integral converge. 

The function $ 3 (xi ; fi) is called the leading- twist distribution amplitude 
(DA). It can be studied using the operator product expansion 7 

^ 3 (x t ; ^J.) = l20f N (fi)x 1 x 2 x 3 ^l + c w {fJ,)(xi- 2x 2 + x 3 ) + Cu(fi)(xi -x 3 ) 

+ c 2 o(m) [l + 7(x 2 - 2xix 3 - 2x1)] + c 2i(M) (1 - 4a; 2 ) (xi-x 3 ) 
+c 22 (m) [3 - 9x 2 + 8x 2 - 12x lX3 ] + ...} (I) 

where (wave function at the origin) and Cik(fJ-) (shape parameters) 

are scale-dependent coefficients which can be defined as matrix elements of 
(multiplicatively renormalizable) local operators. They can be calculated 
using QCD sum rules 8 or lattice QCD. 9 The DA & 3 (xi;fj,) is thus a much 
simpler object compared to the full light-cone wave function tp e ' =a (xi, ki); 
unfortunately this reduction does not seem to work in the Q 2 = 5 — 15 GeV 2 
range. 

Another problem is that the standard pQCD factorization approach 
does not take into account contributions of states with non- vanishing orbital 



a Here and below we do not show contributions which vanish in the limit of small trans- 
verse separations, cf. Rcf. [10]. 
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angular momentum. For example 10 



, - + #^2" 

x | \u t (xi,k 1 )u l (x 2 , k 2 )d l (x 3 , fc 3 )) - \d t (x 1 , k 1 )u l (x 2l k 2 )u l (x 3 , fc 3 ))| 

where k^ = k x ± ik y . The new light-cone wave functions ^{"^(xi, ki) and 
■</4* -1 ( x i; fci) are reduced to next-to-leading twist-4 nucleon DAs 12,13 



$ 4 (x 4 ;m) 



[d?k] 
ttinxz 



J rri]\fX 2 

and, again, can be studied using OPE 7,12 

§i{xi\iJt) = 12Xi(/j,)xiX 2 + l2f N (n)xxx 2 



1+3(1-5x3) 



l + -(l-5a; 2 ) 



+ ■ 



(2) 



Note that to this accuracy twist-4 DAs include one new parameter only, Ai. 
A similar expansion and the reduction to DAs can be worked out for nu- 
cleon resonances. Electroproduction of nucleon resonances at high Q 2 gives 
access to three-quark wave functions, more precisely to the overlap inte- 
grals between the wave functions of the nucleon and the resonance. These 
overlap integrals are in general too complicated to be calculated in QCD 
directly and our strategy will be to reduce these integrals to convolutions 
of distribution amplitudes which can be constrained using lattice QCD. 

We believe that one important physics goal for such studies will be to 
compare DAs (alias light-cone wave functions at small transverse separa- 
tions) of baryon states with opposite parity, e.g. J p = 1/2+ and J p = l/2~. 
It is well known that such "parity-doublets" are non-degenerate in QCD 
because of spontaneous breaking of chiral symmetry. It is not known, how- 
ever, whether this difference mainly affects the "pion cloud" or it is present 
at short distances already and affects momentum fraction distributions of 
valence quarks. Moreover, there are indications 14 that chiral symmetry is 
effectively restored in the spectra of higher-mass resonances and it would 
be extremely interesting to compare the corresponding wave functions. 

The general strategy of combining the constraints on DAs from a lattice 
calculation with LCSRs to calculate the form factors is suggested in Rcf. [15] 
where, as the first demonstration, we considered the electroproduction of 
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Fig. 1. Wave functions at the origin (left panel) and shape parameters of the nucleon 
and TV* (1535) distribution amplitudes (right panel). Only statistical errors are shown. 



iV*(1535), the parity partner of the nucleon. In what follows I describe an 
ongoing work in this direction by the Regensburg theory group and QCDSF 
collaboration. 



2. Nucleon and TV* (1535) Distribution Amplitudes from 
Lattice QCD 

Baryon states of different parity can be identified in a lattice calculation as 
those propagating forward or backward in (imaginary) time, 16 ' 17 so in fact 
the results for iV*(1535) reported in Ref. [15] are essentially a byproduct 
of our calculation of the nucleon DAs. 9 These results were obtained using 
QCDSF/DIK gauge configurations with two flavors of clover fermions for 
two different (3 values and several quark masses on 24 3 x 48 lattices. The 
calculation was done using nonperturbatively renormalized three-quark op- 
erators 18 with up to two derivatives imposing a RF-MOM-likc rcnormal- 
ization condition and converting the results to the MS scheme. In this way 
the mixing with "total derivatives" is automatically taken into account. 

This work will be continued, using larger lattices and smaller pion 
masses in order to minimize effects of the chiral extrapolation. In Fig. 1 
I present preliminary results 19 of the new calculation using ca. 600 Nf = 2 
gauge configurations on a 32 3 x 64 lattice with (3 — 5.29 (a = 0.0753 fm) 
and pion mass ~ 282 MeV (m^L — 3.44)). The left panel shows the 
comparison of nucleon and TV* (1535) wave functions at the origin, /jv and 
Jn* , as a function of mj. The left-most points at ~ 0.08 GeV 2 are new. 
They are much closer to the physical point and show, for the first time, that 
for small pion masses /at* becomes smaller than /jy. This phenomenon still 
has to be understood. 
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Fig. 2. Barycentric plots for the nucleon (left) and N* (1535) (right) distribution am- 
plitudes. 

The results for shape parameters of the nucleon distribution ampli- 
tudes (1) are shown on the right panel. One sees that coefficients of first- 
order polynomials cio and en can be quantified, whereas contributions 
of second order polynomials are less constrained. The main reason for 
this are the 0(a) discretization errors in the chain rule for derivatives 
D(A ■ B) = (DA) ■ B + A ■ (DB) + 0(a) which spoil energy conserva- 
tion: after the (nonperturbative) renormalization we obtain for the sum of 
the quark momentum fractions x\ + xi + X3 ~ 0.94 instead of unity. This 
is one of the issues that have to be addressed in future studies. 

The main result so far is that the accuracy of modern lattice calculations 
is sufficient to detect differences in quark momentum fraction distributions 
in the nucleon and its parity partner state. Our calculations support a 
qualitative picture suggested by QCD sum rules 8 that the valence quark 
with the spin parallel to that of the nucleon carries most of its momentum, 
and for TV* (f 535) this effect appears to be even stronger. As an illustration 
we show in Fig. 2 the DAs of the nucleon and TV* (1535) in barycentric 
coordinates. b Note that the iV*(1535) DA appears to be more narrow and 
shifted towards the lower-right corner. 

All these results are preliminary and the study will be continued using 
larger lattices and smaller pion masses. The comparison of the nucleon and 
iV*(1535) will remain our primary goal for some time, but the calculations 
will also be extended to the full J p = l/2 + and J p = l/2~ baryon octets 
and to the decuplets. The main problem that has to be addressed in fu- 



b In difference to the similar plot in Ref. [15] we only include the terms in cio and cxi 
and discard contributions of second order polynomials which have large errors. 
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Fig. 3. Schematic structure of the light-cone sum rule for baryon form factors. 

ture is the identification of resonance contributions for small quark (pion) 
masses such that strong decays, e.g. TV* (1535) — > Nn,Ni], are allowed. 
The separation of TV* (1535) and iV*(1650) may prove to be difficult. In 
particular the large decay width of iV*(1535) in the Nr] channel has to be 
understood and it can only be addressed in lattice calculations with three 
flavors of dynamic fermions. Discretization errors 0(a) in the lattice defi- 
nition of the relevant operators become a serious issue for second moments 
and have to be reduced in order that the results for C2k coefficients become 
fully quantitative. 

3. Transition form factors from Light-Cone Sum Rules 

The matrix element of the electromagnetic current j™ between spin- 1/2 
states of opposite parity can be parametrized in terms of two independent 
form factors, which can be chosen as 



where q = P' — P is the momentum transfer. The LCSRs are derived from 
the correlation function 



where r\ is a suitable operator with nucleon quantum numbers, e.g. the Ioffe 
current. 20 Making use of the duality of QCD quark-gluon and hadronic de- 
grees of freedom through dispersion relations one can write a representation 
for the form factors appearing in (3) in terms of the DAs of N* . Schemati- 
cally, the sum rules take the form 



(N*(P')\f™\N(P)) = u N ,{P')^T v u N {P) : 




(3) 




Gia{Q 2 ) = ^ / [dx] C 1}2 ;k(xi, Q 2 , s , M 2 , fi, a s (fi)) $k(%i, m) 
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where the sum goes over contributions of nucleon DAs of increasing twist 
and Ck are the coefficient functions that can be calculated in QCD per- 
turbation theory and modified using dispersion relations to include two 
nonperturbative parameters: So, the interval of duality, and M 2 , the Borel 
parameter which specifies the distance (in imaginary time) on which match- 
ing between hadronic and quark representations for the correlation function 
is being done. The dependence on M 2 is unphysical in the same sense as 
the factorization scale \x dependence in truncated perturbative expansions, 
but it is usually weak. The pQCD limit 3-6 corresponds to the leading part 
of the contribution of leading-twist DAs at large momentum transfer, and 
the main difference is that higher-twist contributions are not suppressed by 
powers of Aq CD /Q 2 , but rather of Aq CD /s with s ~ (1.5 GeV) 2 . The 
reason for this is that LCSRs take into account "soft" contributions to 
the form factors coming from large transverse separations. The attractive 
feature of this approach is that such terms are calculated in terms of the 
DAs, thus avoiding the need to know (model) the full nonperturbative k± 
dependence of wave functions. 

In leading order, the sum rules for Q 2 G\{Q 2 ) / '(toattojv*) and —2G 2 (Q 2 ) 
have the same functional form as the similar sum rules 21,22 for the Dirac 
and Pauli electromagnetic form factors of the proton, with the replacement 
Tti j\[ — > tun* in the light-cone expansion part, and different DAs. 

The experimental results for the electroproduction of A* (1535) are usu- 
ally presented for helicity amplitudes Ai/ 2 (Q 2 ) and Si/ 2 {Q 2 ) which can be 
expressed in terms of the form factors: 23 



Here e is the elementary charge and B, C are kinematic factors defined as 



The results of the LCSR calculation of the form factors (normalized to 
the dipole) and helicity amplitudes using lattice-constrained N* DAs from 
Ref. [15] is presented in Fig. 4. The shaded areas show the estimated un- 
certainties. The Q 2 -dependence of the form factors Q 2 G\{Q 2 ) and G 2 (Q 2 ) 
is predicted to be similar to the Dirac and Pauli nucleon electromagnetic 
form factors, respectively. Only the normalization is different: The observed 
negative amplitude Si/ 2 implies a small value of G 2 . 



A 1/2 =eB Q 2 G 1 (Q 2 ) + m N (m N , - m N )G 2 (Q 2 ) , 
S 1/2 = -^=B C \{m N ~ niN^G^Q 2 ) + m N G 2 {Q 2 ) . 
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Fig. 4. Form factors (left panel) and helicity amplitudes (right panel) for the electroexci- 
tation of the AT* (1535) resonance calculated using the DAs in Ref. [15]. The experimental 
data are from Refs. [24-27]. 



These results have to be viewed as exploratory and further work is 
needed to make them quantitative. First of all, the LCSRs for baryon form 
factors have to be extended to the next-to-leading (NLO) accuracy, i.e. 
including the 0(a s ) corrections. This calculation is rather cumbersome and 
not straightforward because of contributions of evanescent operators. As the 
first step in this direction, the NLO corrections to contributions of leading- 
twist DAs in the LCSRs for the electromagnetic nucleon form factors Fi 
and F2 have been calculated in Ref. [28], see Fig. 5. One sees that such 
corrections are significant, especially in the Ge/Gm ratio. The extension of 
this calculation to contributions of sub-leading twist-4 DAs and the detailed 
study of baryon mass corrections to the sum rules is planned for the coming 
years. 

Second, one has to remember that the LCSRs suffer from irreducible 
uncertainties due to the duality assumption for contributions of higher res- 
onances and the continuum which for meson form factors is estimated to 
be of order 10%. (see e.g. Ref. [29]). The experience with applying this 
method to baryons is much less than for mesons, so that this "systematic 
error" cannot be estimated reliably. It is therefore imperative to apply the 
same technique to a maximally broad class of reactions. At present applica- 
tions include nucleon electromagnetic and axial form factors, 21 ' 22 ' 30 NjA 
transitions, 33 pseudoscalar- and vector-meson couplings to octet and decu- 
plet baryons, 34 weak decays of the type A& — » pivg 31 and Aj, — » A£ + £~ , 32 
and pion electroproduction at threshold. 35,36 I expect that this list will 
continue to grow, giving us confidence in the whole program. 
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Fig. 5. LCSR results for the electromagnetic proton form factors for a realistic 
model of nucleon distribution amplitudes. Left panel: Leading order (LO); right 
panel: next-to-leading order (NLO) for twist-three contributions. Figure adapted 
from Ref. [28]. 
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